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Abstract 



i-C I We study the exclusive decays of -B K^*H^l~ within the framework of the PQCD. 

We obtain the form factors for the B K^*^ transitions in all allowed values of g^, which 
agree with the lattice results. We find that our distributions of the decay rates and leptonic 

^ . asymmetries are consistent with that given in the other QCD models in the literature. 
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1 Introduction 



The recent CLEO measurement of the radiative b ^ s'j decay has motivated theorists to 
study exclusive rare B meson decays such as B ^ K^*H~^l~ 0. In the standard model, these 
rare decays occur at loop level and provide us with information on the parameters of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements as well as various hadronic form 
factors. In this paper, we examine the decays of B ^ {K, K*)l^l^ within the framework of 
the perturbative QCD (PQCD). 

The calculations of matrix elements for exclusive hadron decays can be performed in the 
PQCD approach developed by Brodsky-Lepage (BL) The application to B meson decays 
was first carried out in Refs. |^ and [Q. In the BL formalism, the nonperturbative part is 
expressed as the hardon wave functions which could be determined via various QCD models 
such as the QCD sum rule method or lattice gauge theory and the transition amplitude 
is factorized into the convolution of hadron wave functions and the hard amplitude of the 
constituent quarks. However, with the BL approach, the nonperturbative effects appear 
1^, 1^ if one of the constituent quarks carries nearly all the momentum of hadron. To 
solve the problem, Li and Sterman proposed by including the transverse momentum of 
constituent quark, kx, and the Sudakov form factor to the wave functions to suppress the 
soft contributions from higher order corrections. In terms of the parameter b with b being 
the conjugate variable of kx, they also showed that the effects can be also expressed as that 
in the BL factorization formalism. 

The modified PQCD factorization theorem for exclusive heavy meson decays has been 
developed some time ago ||TU|, [TT|, and applied to nonleptonic B D^*\{p) ||T^, penguin 
induced radiation B K*'~f [|l3l, and B KK decays. These decays involve three 
scales: the My/ scale as the initial condition of renormalization-group (RG) equation, the 
typical scale t which reflects the specific dynamics of the heavy meson decays, and the fac- 
torization scale 1/6. Above the factorization scale, there are two large logarithms \\i{My/ /t) 
and ln(t6), generated from radiative corrections. The former gives the evolution from Mw 
down to t described by the Wilson coefficient (WC), while the latter from t to 1/6. There 
also exist double logarithms ln^(P6) arising from the radiative correction to the meson wave 
function, where P is the dominant light-cone component of a meson momentum. Resuming 
these double logarithms leads to a Sudakov form factor of exp(— s(P, 6)) which suppresses 
the long-distance contributions in the large b region such that the applicability of the PQCD 
around the energy scale of the bottom quark mass could be guaranteed. 

The typical three-scale factorization formula is generally written as the following convo- 
lution product 



C{t)®H {t) ® (p (x, b) ® exp [s (P, 6) - 2 /* ^7 (a, [jj))] (1) 

V h/b fx J 



where C(t), H{t) and (f){x,b) denote the WC, hard decay amplitude and non-perturbative 
wave function, respectively, and the quark anomalous dimension 7 = — ag/vr is evaluated from 
t to 1/6. Except (t){x,b) dictated by non-perturbative dynamics, all the convolution factors 
in Eq. (|l]) are calculable. Note that differing from the conventional factorization assumption 
(FA), the WC is also one of the convolution parts in Eq. (|l|). Thus, the /i dependent problem 
occurring in the FA could be solved naturally in the three-scale factorization formula. 

The paper is organized as follows. In Sec. 2, we study the form factors in the framework 
of the PQCD for the decays of P — ^ K^*'> transitions. In Sec. 3, we derive the forms of 
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the differential decay rates and lepton asymmetries for B K^*H'^l^ based on tlie PQCD. 
In Sec. 4, we give tlie numerical analysis. We will also compare our results in the PQCD 
approach with that in the other QCD models. In Sec. 5, we present our conclusions. 



2 Form factors in the framework of the PQCD 



In the decay of i? ^ Hl^l , the momentum of B [H) in the light-cone coordinate is chosen 
as Pi = (Pi"J:2)'^i(2)'0t), where Pf = Mb/V^ and = {Eh ± Ph)/V2 with Eh = 

(M^ + Mfj — q^) /2Ms, Ph = \/ ^h — Phi ^"^^ 'f is the squared momentum transfer. We 
define the momentum of the light valence quark in the B meson as ki and use xi = kf / Pi 
with ki and kiT being the plus and transverse components of /ci, respectively. The two light 
valence quarks in the H meson carry the longitudinal momenta X2P2 and (1 — X2)P2, and 
transverse momenta \i.2T and —\<.2t, respectively. 

The Sudakov resummations of large logarithmic corrections lead to the exponential forms 
of exp(— 5*5) and exp^—Sn) for B and H wave functions, respectively, where 



Ssit) = .(xiP+,6i) + 2 r ^7(«s(/x)) 

Jl/bi jjL 



Suit) 



s{x2P^,b2) + s{{l-X2)P^,b2) 



2f %{aM). 

Jl/b2 jJ 



(2) 



In Eq. (|^), 6iand 62 represent the transverse momentum extents of B and H and are 
conjugate to the parton transverse momenta kiT and k2T, respectively. The form for s is 
written as [^, |T^, |^ 



(3) 



where the anomalous dimensions A and D calculated to the two and one-loop levels, respec- 
tively, are given by 



A 



D 



TT 



2«, 

3 TT 



In 



67 _ tH 

,27i5-l\ 



27'^ 3 V 2 



TT 



(4) 



with Cf = 4/3 being a color factor, / = 4 the active flavor number, and 7^; the Euler 
constant. For the running coupling constant, we use 



47r 



/3oln(/iVA|cD) 



(5) 



with /?o = (33 - 2/)/3. 

To get the transition elements of B ^ H {H = K, K*) with various types of vertices, 
we parametrize them in terms of the relevant form factors as follows: 



(ir(P2)| VjP(Pi)) = Pi g 



P^. + F2[q- 
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with 



Friq 



{K*{P2,e)\V,\B{P,)) = e^_^^-pV, 
{K*{P2,e)\A^\B{P,)) = Ao{q')e; + e*-q[A,[q')P^ + A2{q')q, 



{K*{P2,e)\T^,q'' \B{P,)) = (g^) e^,„^£-pV, 
K*{P2,e)\T;:X \BiPi)) = -To{q')s;-e*-q[T^{q')P, + T2{q' 



To(g2) + [Ti(g2)p.g + T2(g2) q' 







(6) 



(7) 



s ia^yh, and T^^ = s ia^y'-)^ b. The correspondences 



where = sj^b, A^ = s 7^75 b, T^^ 
of our notation to that usually used in the literature are shown in the Appendix. Using the 
PQCD factorization formula, the components of form factors defined in Eq. (P) are found 
to be 



Fi(g2) = -MCpMl j^[dx] p bidbib2db2 (t)B{xiM) 
{[r'K (2a2i^ + 2P2K - 1) <P'k {x2, h) 

- (1 + f32K + (1 + 2rK - s) a2K) (pK {X2, fe2)] Ek {i^^) {xi,X2, 61, 62) 
+ [-2r^ (1 - aiK - Pik) <\>k (2^2, ^2) 

+ {r-K (1 - OL^K + Pik) - s Pik) <Pk (^2, &2)] Ek {fT) h'' {x2, ^i, 62, &i)} (8) 

F2((f) = -SttCfMI Ijdx] 1^ bidbib2db2(pB{Xl,bi) 

{[-r'Kil + 2a2K - 2P2K)<P'Kix2, b2) + (1 + P2K 

+ (1 + 2rK - s)a2K)(l)Kix2, b2)]EK (ti^^) (xi, X2, 61, 62) 

+ [2r^ (1 - aiii- + Pik) 0'^ (x2, ^2) + {rKi^ - ^iK - Pik) 

-(2 - s)P^K)(t^K {X2. h2)] Ek (tf ) h"" (x2, Xi, 62, &i)} (9) 

FT(g^) = -SuCfMb J^[dx] bidbib2db2(pB{xi,bi) 

X { [-a2x?^^0ii- (a;2,&2) + (1 - 2P2k) 4>k {x2M)\ Ek (4^^) {xi,X2MM) 
+ [2r;, (1 - aiK) <P'k (3^2, &2) - PiK(t>K {X2, h2)] Ek (tf ) h"" (x2, Xi, 62, &i)}(10) 

K (g^) = SuCfMb [dx] bidbib2db2 0b (xi, 61) (j)K* (^2, ^2) 

X {[1 - 2/32/r* - ^*a2K'] Ek' {t^P) h'^* (xi, X2, 61, 62) 
+ [v^(l-«ii^*)]^i^' (ti'))/i^*(x2,xi,62,&i)}, (11) 
^0 (?^) = SuCfMb [dx] bidbib2db24>B {xi,bi) (pK* {x2,b2) 

X {[(1 + Tk* - s) (1 - 2/52i^. + y/r]^a2K*) 

-ArK'a2K'+2^,{l + P2K')]EK* {t^^'>) h""* (xi, X2, &i, 62) 

[(1 + ^i^* - ^) (1 - «iK*) - 2Ai^*] ^i^* (^i'O ^^2, xi, 62, &i)} (12) 
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= SttCfMb [dx] bidbib2db2 (f>B{xi,bi)(f)K*{x2,b2) 
X {[-1 + 2P2K* + Vr^0(2K*] Ek* (ti^^) h^" (xi, X2MM) 
+ [v^(-l + 2P^k*)]Ek* [tf^) h""* {x2.x^MM)] , (13) 

^2 (?^) = SttCfMb [dx] bidbib2db2 4>b {xi, bi) 4>k* {x2, ^2) 

X {[1 - 2P2K* - V^a2K*] Ek^ h""* (Xi, X2, 61, 62) 

+ [v^(l - + 2/9ii^0] (^i'O (^2, x^, b2, 61)} , (14) 

T = SttCf M2 ^'[dx] bidbib2db2 (pB {Xl, 61) (X2, 62) 

X {[- (1 + V^) + 2sQ;2i<:* + (2^^ - 1) (q;2k* + P2K*)] 

XEk* (4'^) /l"^* (xi,X2,6i,62) 

["IK' + Ai^* (xi) - (42)) /i^^* (a;2,a;i, 62,^1)}, (15) 

To (g') = -SttC^M^ ^'[dx] bidbib2db2 (t>B (Xl, 61) (X2, &2) 

X {[(1 - tk* - s) (1 + /32i^*) + (1 - ^K* + s) 

-2v/r^((l - tk*) {a2K* + /32ii:*) + s {I32k* - a2K*)) 

-2rK*a2K- + (1 + rK* - s) (1 - s) a2K-] Ek* (4^^) (2^1, X2, &i, 62) 

[(1 - OiiK* - /3lK*) (1 - Tk*) 
-S (1 + PiK* - OliK*)] Ek* (t^e^^ (X2, Xl, 62, bi 

)}, (16) 

(?^) = -StiCfM'^J^ dxidx2 bidbib2db2 (f)B {xi,bi) (pK* {x2,b2) 

X {[.sa2A'* - (1 + y^) - (1 - 2Vrf?^) («2K* + I32K*)] 
xEk* (2:i,2:2,&l,&2) 

[1 - aiK* - PiK*] Ek* (^i'^) h""* {x2, x^, 62, 61)} , (17) 
^2(5^) = -SnCpM'^ dxidx2 &i<i&i&2'^&2 0b (a^i, &i) (2^2, ^2) 

X {[{2rK* - S) a2K* + (1 + Ci2K* + /92i^-) 

- (1 + 2a;2X* - 2(52K*)] Ek* {t'-p) /i^* (xi, X2, 61, 62) 

[1 - aiK* + M Ek* (if) /i"^* (^2, x^M. &i)} , (18) 

where 0s, 0/^ (0'/^)) and (j)K* are the wave functions of B, pseudovector (pseudoscalar) of 
K, and K* mesons, respectively, the evolution factor are given by 

Eh (t) = «s (t) exp {-Sb (t) - Sh (t)) , (19) 
and the related kinematic variables are parametrized as 

aiH = ^Xi, (j^H = -il^ — Xl, 

2 \ J 

a ^ ( I ^^^^H - s\ 

P2H = 1=X2, Ot2H = o IH ^ X2, 

JTp^i 2 \ / 



M, 



B 



Ml' 



(20) 



with 



'l-rHf-2s{l + rH) + s\ 
Ml 



rriQK -- 

nis + rrid 

The hard functions, , are written as 

{xi,X2MM) = Kq {DHy/xiX2bi) [9 {hi - 62) Ko {Dn^/x^hi) Iq (-Dj/v^^a) 



(21) 



f22) 



with 



D 



Ml 



H 



Ml 



The derivation of /i, from the Fourier transformation of the lowest-order hard decay amph- 
tude, is similar to that for B KK decays [|I^. The hard scales t*^^'^-* are chosen by 



max {y/x^Dn, l/h, I/62) 
max (y^Dn, I/&1, 1/&2) 



(23) 



The wave functions 0^ and 0^ are defined by |jT 



P. 



— 'TK 



X 



dy+ 

271 

dy+ 

277 



\u{y^h-i,s{o)\H) 



-ixPo y 



(24) 



with the normalization conditions of 



dx(j)B{x) 



dxcpH^x) 



dx(f)'j^{x) 



fB{H) 



JO JO 2^2Nc 

We note that unlike the kaon case, we do not distinguish the pseudovector and pseudoscalar 
components of the B wave functions since the factor MB/ijrib + m^) is close to one. We also 
note that from Eqs. (|T5|) and (p!7| ) we obtain T(0) = — Ti(0). 



3 Differential decay rates and lepton asymmetries 

The effective Hamiltonian of 6 ^ s /~ is given by pO| 



GpaXt 



TT 



(25) 



where Cs{fi), Cg and Ct{h) are the WCs and their expressions can be found in Ref. |^ 
for the SM. Since the operator associated with Cg is not renormalized under the QCD, it is 
the only one with the /i scale free. Besides the short-distance (SD) contributions, the main 
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effect on the branching ratio comes from cc resonant states such as ... etc., i.e., the 

long-distance (LD) contributions. In the hterature |22|, |2^, it has been suggested 
by combining FA and vector meson dominance (VMD) approximation to estimate LD effects 
for the B decays. Hence, including the resonant effect (RE) and absorbing it to the related 
WC, we obtain the effective WC of Cg as 



3 _ , vrr (j ^ /+/-) M, 

Kj 



(26) 



where we have neglected the small Wilson coefficients, h{x, s) describes the one-loop matrix 
elements of operators Oi = Sa^^Pi^ii cpJ^PiCa and O2 = s'j'^PLb c'j^Plc Mj (Tj) are 
the masses (widths) of intermediate states, and the factors kj are phenomenological param- 
eters for compensating the approximations of FA and VMD and reproducing the correct 
branching ratios Br {B J/^X l+l-X) = Br {B ^ J/^X) x Br (J/^ l+l'). In 
this paper, for simphcity, we take kj = —1/ (3Ci (/i) + C2 (/x)). 

Using Eqs. (^) and (p5D, the transition amplitudes oi B —>■ {K, K*) l~ are as follows 



(27) 



and 



M 



GpaXt 



K* 



2v/27r 



(28) 



In Eqs. 



iV\q')e,,^pe*''P^q^ - Al{q')e; - e* ■ qAl{q')P,] lYj.l} . 
and (pSD, we have included the WCs by inserting them into Eq. ( [TI3| ) as 

Ei (t) = C, (/i) as (t) exp {-Sb (t) - Sh (t)) (29) 



and the superscripts of form factors denote the associated WCs of Cg^'^ , Cg, and Cj, with 
the new definition of Ejj in Eqs. (||)-(|1^), respectively. 

As usual, after integrating the angle dependent phase space, the differential decay rates 
for B Hl+l- {H = K, K*) are found to be 



dTn {q^ 
ds 



G],a^\\t\"Ml 
3 X 2^5 



\ 



where 



F^ (g2) - 2m,Fl [q' 



Amf 



f! (g^) 



1 + 



2vai 



mf 



13. + 12^* 



H 



Ft [q') f + (1 - Tk) Re F^ (g^) Ff (g^) + |f| (g 



s {2^K*V (g2) + 3Fo (g^)) + (f, (f) + ^K*h (f) 



+2 (1 - TK* - s) Foi (q 



(30) 

(31) 
(32) 

(33) 
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-2 



V (g^) Ms 



Mf 



+ 



Al {q') Ms 



2 {l+TK*) - S 



1 - 

2r 



A'* 



Re A? (g')MBA^* (g^) 



(34) 



with 



l/s - ^!^T^ (g^) +\v\q-)Ms 



q 

2mbT^{q^) 



B 



q' 



+ 



Al iq'' 



At {q') 
Re 



2mbMs 



Ms 

2 



q 



2 TlU 



Al (g^) Ms 



+ Re 



Ms q- 
Al {q') 



Mr. ^ J ' 



IB 



LB 



(35) 



The forward-backward asymmetry (FBA) can be defined by 

1 



A 



FB 



dV (s) /ds 



1 . d^Vis) 
d cosu 







dsd cos 9 



C^COS^ n 

-1 dsd COS tl 



(36) 



where 6 is the angle of charged with respect to the B meson in the rest frame of the lepton 
pair. For B K*l^l^ decay, the FBA is found to be 



•^FB 



VA 



1 + '-^]Pk'+ Vl^^K* 



(37) 



with 



Re 
+ Re 



M, 



B 



V' [e) Ms - '-^r (g^) 



AW) 2m,Tj(g2)\ . . 

-w^ T^^irr i^J^^ 



B 



q^ 



M 



B 



(38) 



As expected, the FBA in Eq. (|37| ) is sensitive to the chiral structure of interactions since 
it is related to the product of V and A currents. It is clear that the FBA for B — > Kl^l~ 
vanishes since there is no form factor from the axial current. 

Another interesting lepton asymmetry is the longitudinal polarization of the lepton, de- 
fined by 



dr(n=-l) _ dr(n=l) 
ds ds 

dr(n=-l) _|_ dT(n=l) 
ds ds 
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where n is the projection of the lepton / momentum to the spin direction in its rest frame. 
For B {K, K*) l^l", the polarization asymmetries can be expressed as 



and 



with 



Re 



F^{q')-2m,F^{q'))Fr {q' 



2^/l-^ 



{s [2^K*Rv (g') + 3i?Ao {q' 



4r 



Rao (q^) + Vk*Ra^ (q^) + (1 - r-A- - s) Ra,,, [q' 



(39) 



(40) 



Rv 



^^0 W 



RAi 
Raoi 



Re 
Re 
Re 
Re 



y8 (^2) - ^!^T^ (g^) ) (g^) Ms 



q 



, Mb 

(g^) M, 



g2 Mb J Mb 
2mbMB 



Tl (g^) At (g^) 



'AW) 2m,Tj(g2)\ . . 



+ [Al{q^)MB-^-^Tl(q-) 



q^ 



Mf 



respectively. 



(41) 



4 Numerical Analysis 



4.1 Form factors 

In Eq. (|l]), 4>{x,b) is the universal wave function and cannot be calculated perturbatively. 
However, due to the universality, we can determine it by matching with the B decay exper- 
imental data. With the ratio of 



R 



B{B^ KH^) 



0.95 ±0.3 



(42) 



given by the CLEO measurement |^, where B{B^ K^tt^) represents the CP average of 
the branching ratios B{B^ — * K^tt~) and B{B^ K^tt^), one can get the proper wave 
functions <^b, (px, and (f)'j^ ||T^ while (f)K* can be done by the branching ratio of i? ^ K*'y 
||r3| . For the B meson wave function, we take 

1 /xMf ^2 , .2 U2- 



')b{x^ b) = Nbx'^{1 — xY exp 



Ub 



(43) 
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with the shape parameter ujb = 0.4 GeV [^. The normahzation constant Nb = 91.7835 
GeV is related to the decay constant fs = 190 MeV. The kaon wave functions are chosen as 

(t)K{x) = -^/i,x(l -x)[l + 0.51(1 -2x)+ 0.3(5(1 - 2x)2 - 1)] , 
3 

Mx) = -^/i,.x(l -x)[l + 0.51(1 -2x) + (5(1 -2x)2-l)], (44) 



where (^k is derived from QCD sum rules 2a], and the second term in the expression of 



(j)K corresponds to SU{3) symmetry breaking effect. The decay constants fx and fx* are 
set to be 160 and 190 MeV (in the convention of /^r = 130 MeV), respectively. Note that 
the intrinsic b dependences of wave functions in Eq. (0) are neglected. However, this is 
a good approximation only for the fast recoiling meson, that is, the transverse extent of 
the wave function is less important in the energetic outgoing situation. With above wave 
functions and taking Mb = 5.28 GeV, Mk = 0.49 GeV, m, = 100 MeV and Mk* = 0.89 
GeV, the form factors of B K defined in Eqs. (P) as a function of are shown in Figure 
1. The values of the form factors at = are given in Table 1. We now compare our 
results with that in the light cone-QCD sum rule (LCSR ) [Q. Using the identities in the 
Appendix, we find that except V{0) is slightly smaller than that of the minimal value, while 
the remaining form factors are within the allowed values, in the LCSR. Recently, it has been 

Table 1: Form factors at = in the PQCD. 



Fi(0) 


^2(0) 


Ft(0) 


ViO) 


^(0) 


^i(O) 


^2(0) 


T(0) 


Ti(0) 


7^2(0) 


0.33 


-0.267 


-0.054 


0.063 


2.02 


-0.05 


0.059 


-0.350 


0.350 


-0.281 



mentioned that by combining large energy effective theory (LEET) [T^ , originally proposed 
by |TB[, with the measurement of 5 — > K*j, the form factors V^(0) and Ao{0) could be fitted 



model-independently to be [32 



V{0) ~ 0.069 ±0.011, 

Ao(0) ~ 1.650 ±0.114. (45) 

Hence, from Table 1, we clearly see that the values of V{0) and Aq{0) are within la and 
3.2a of values in Eq. (|45|) , respectively. It is worth to mention that the LEET predicts 
T2(0)/y4i(0) ~ 4.89 while that in our approach is 5.62. For the comparison of the form 
factors between different models at = 0, one can refer to Ref. |^ for a more detail 
analysis. 

For the exclusive B — > K^*H~^l~ decays, if we use b independent wave functions, as the 
mesons reach the slow recoil, the suppression in the large b region is weaker such that form 
factors will blow up at points away from = as seen from Figure 1. It is inevitable to 
include the intrinsic b dependence to the outgoing meson wave functions. However, the effect 
of such dependence is less significant for the small g^ values than that of large ones. Instead 
of using an exponential b dependent form as the one for the B meson wave function in Eq. 
(H3), for simplicity, we use the trial wave functions as 
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0S?(x,g2) = (1 - ^)^f+^/ll (^) • (46) 

On the other hand, since the available region of the PQCD essentially cannot include all 
allowed values of g^, to have the form factors in the whole accessible values of g^, we would 
adopt the parametrization of effective form factors as follows: 

F{S)=F (0) exp ((TiS + (725^ + (TgS^) (47) 

with s = /M"^ to fit the values up to ~ 15 GeV^ calculated by the PQCD. By extrap- 
olating to near the endpoint of g^, we found that the values of form factors are consistent 
with lattice results |2^, 0. To illustrate how good the trial functions in Eq. (|46| ) are, we 
show the form factors for i? ^ i^' in Figure 2. From the figure, we find that our results are 
basically the same as that from the QM and LCSR 0. 



4.2 Decay rates 

With the confidence of calculating the form factors by using Eqs. (|46| ) and (0), we now 
study the decay rates oi B K^*U^l~ . Unlike the conventional FA, the WCs in Eqs. (|27| ) 



and (^) are the members of integrations in the PQCD. Thus, adopting the approach similar 
to form factors, we calculate the transition amplitudes with the wave functions in Eq. (|46| ) 
and the exponential forms in Eq. (^) to fit the values calculated by the PQCD for the whole 
range of g^. After integrating g^ dependence in Eq. (^), the decay branching ratios without 
including LD contributions for B —>■ {K, K*) are listed in Table ^ and their distributions 
for the differential decay rates are shown in Figure 3. Comparing with the curves in the QM 
and LCSR, we find that the differential rates of i? — > Kl^l~ are consistent with each other. 
However, there exists a slight difference in S ^ K*l'^l~ . There are two main reasons for the 
difference: (a) /i scale dependence for the WC and (b) the effects from Ai(g^) and A2(g^) in 
the PQCD. For the results in QM and LCSR, we have used the WC at /i ~ m;, as done in the 
literature, whereas that in the PQCD, /i scale is the typical scale t determined by Eq. (p3[) . 
As seen from Eq. (^) the effect of ^i(g^) is large since there is a factor of Mb associated 
with it, while that of ^2(9^) only affects in the mode of -B ^ K*t^t^ due to the lepton 
mass dependence. Thus, measuring the exclusive modes of B K*l^l^ would distinguish 
various QCD models due to the difference shown in Figure 3. From Table ^ we see that our 
PQCD results of the decay branching ratios for B — > K*fi~^fi~ and B K*e^e~ are quite 
different. This can be understood by noting that in Eq. (Q) there is pole of g^ associated 
with the photon penguin induced couplings. These pole terms make the rates sensitive to 
the kinematical region of g^ > 4m?. 



4.3 Forward-backward asymmetry 

^From Eq. (P7|), we present the forward-backward dilepton asymmetries of -B ^ K*l^l^ [l = 
fi, t) in Figure 4. We note that the FBA of B ^ K*e~^e~ is similar to that of the muon 
mode in Figure 4. For the light lepton pair decays such as B K*fi~^fi~, we see that the 
FBA is positive at low g^, gets zero at q^ /M"^ ~ 0.16, and then becomes negative. Due to 
the large terms related to Ai[q^), the rate at lower g^ in the PQCD has a larger value than 
that in the QM and LCSR. As mentioned in Ref. |0] with the FA, the location of zero point 
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Table 2: Deacy branching ratios in the various QCD models without including LD effects. 



Mode 


PQCD QM LCSR 


WB{B KI+1-) 
WB{B Kt+t-) 


5.33 5.56 5.20 
1.29 1.28 1.25 


106B(fi ^ K*e+e-) 
IQ^B{B K*i2+i2-) 
10^B(fi ^ K*T+T-) 


2.26 1.88 2.23 

1.27 1.49 1.78 
1.24 1.43 1.77 



is only sensitive to the WC and insensitive to the form factors. However, since with the 
three-scale factorization formula the WC cannot be factored out of the transition amplitude 
and it is uncertain to choose the universal wave function of K*, the determination of zero 
point is harder as the test of the SM in the approach of the PQCD, unless we can fix K* 
wave function more precisely. On the other hand, it is worth to mention that when the 
sign{C^Cj) = +, opposite to the SM, the zero point disappears. Thus the FBA is quite 
sensitive to the sign of the WC and can be used as a good candidate to test the SM. 

4.4 Polarization asymmetry 

The lepton polarization asymmetries of B —* {K, K*) l+l^ are displayed in Figure 5. For 
B Kl+l~ , Pl is equal zero at = and g^|max = {Mb — because it is related to 

y^l — Amf /q'^ifK- Without LD effects for the light leptons (/ = e, /i), from Eqs. (|3TD and 
(p9|), we easily realize that Pl is near —1 in the most region of g^. However, for B —>■ K*l+l~ , 
since ipx* cannot be factored out in the numerator, there is no vanishing point at g^|max; 
and the transition matrix elements of the set {T'^} are always associated with a pole of 
q^. Hence, at the low momentum transfer region, penguin induced electromagnetic effects 
are dominant. Also comparing Eq. ( ^OD with Eq. (0), the related terms of {T^} are two 
powers for the differential decay rate but only one power for the numerator of Pl so that the 
magnitude of the distribution at low has a smaller value. From Figures 5c and 5d, the 
polarization asymmetry of B ^ K*fi+fi~ in the PQCD has slightly different distribution to 
other models at low g^, especially that the deviation for B — > K*t+t- is large. However, 
according to Figure 6, we find that our results are comparable with that given by the light- 
front formalism (LF) []33| , As mentioned before, the influence of both larger values 
from the Ai (g^) and A2 (g^) terms in our approach is visible for the r mode. Therefore, 
by measuring the longitudinal r polarization in B K*t+t~, we can either determine a 
more proper K* wave function or test the feasibility of our PQCD approach for semileptonic 
decays. 

5 Conclusions 

By the three-scale factorization theorem, we have gotten the form factors in g^ < 15 GeV^; 
and with the parametrization in terms of the exponential forms to extrapolate the g^ depen- 
dent form factors to g^ax; have obtained the consistent results with that from the lattice 
PP| , With the PQCD, we have pointed out that the largest uncertainty in our results 
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is from the non-perturbative wave functions. Though the universal wave function could be 
determined by some non-leptonic decays, the intrinsic b dependence which suppresses the 
soft dynamics contribution is still unknown. 

With the kaon wave functions fixed by the decays of B ^ ttK and assumed depen- 
dences, we have shown that the distributions of the decay rates and leptonic asymmetries for 
B — > Kl^l^ in the PQCD agree well with that in the other QCD models such as the QM and 
LCSR. However, there are some differences for that in S — > K*l'^l~ among the various QCD 
models. Finally, we remark that although B — > K*^ could also give us some information 
of the K* wave function, one still cannot fix it satisfactorily due to various uncertainties 
in the decay. Moreover, the assumption of the same q"^ dependent factors in B ^ K is 
not necessary for B K* and thus, to have a reliable calculation, we need more precision 
measurements involving the vector kaon meson in order to settle down the K* wave function. 
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Appendix 

In order to connect our form factors in Eq. (|^) to those usually used in the literature 
0, ^ in this Appendix we show explicitly the relationships among the form factors. 

In terms of the notation in Refs. 0, |3^, the form factors for B {K,K*) decays with 
respect to various weak currents are parametrized as 



{K{P,)\V,\B{P,)) = U q') P, 



P-q 



q' 



P • q 



q' 



{K{P2)\T,,q' \B{P{)) = -{P,q'-q,P-q)fT{q') 

V iq') 



{K* iP2,e)\V, \BiP,)) 
{K* {P2,e)\A, \B{P,)) 



*u r>a 3 



Mb + Mk' 



(48) 
(49) 

(50) 



+i{MB + MK*)A^(q^] \ e 
e* ■ q 



Pu 



P-q 



{K* iP,,s)\T,,q'' |5(Pi)) 
K* {P2,e)\r^X IS (Pi) 



n (g'j e,.af3e*''P''q 
tT^{q') [e;p.q + e*-qP, 



Mb + Mk' \' ^ q" 

*U TDCl B 



-tT' q')e*-q 



P-q 



(51) 
(52) 

(53) 



where V^= s 7^ 6, = s 7^,75 6, T^.^ = s io^^h, Tj^^ = s ia^^lb 6, P = Pi + P2, g = Pi - P2, 
and P -q = — ML,) . Redefining the wave functions and comparing to Eq. (P) , we obtain 



Pi 
F2 

V 
Ao 
A, 

A2 
T 
To 

Ti 
T2 



u 

Ml-Ml 

2 Uo — J+) 

v 



Mb + Mk* 
{Mb + Mk*)A\ 

A'2 

Mb + Mk* 

4 [2M^*A[, - (Mb + Mk*) + {Mb - Mk*) 



q^ 



T' 



-T' 
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Here we have neglected to show the dependence for the form factors. From the above 
identities, we find some interesting relations at = and they are given by 



/o(0) = /+(0), 
Ti(0) = T^(0). 

From Eqs. (|15D and (p!?!), we get T(0) = — Ti(0). Hence, based on the modified PQCD 
factorization theorem, we obtain the relation T[ (0) = Tl^i^]) that is the same as that in Eq. 
(3.6) of Ref. 0. 
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Figure Captions 

Figure 1: Form factors of Fi (solid curves), F2 (dashed curves) and Ft (dotted curves) 
for the B ^ K transition as a function in the PQCD (bold lines) and QM 
(un-bold lines), respectively. 

Figure 2: Form factors of Fi, F2 and Ft for the B ^ K transition as a function with 
the dependent wave function in Eq. ( ^6|) in the PQCD (solid curves), QM 
(dot-dashed curves), and LCSR (dashed curves), respectively. 

Figure 3: The differential decay braching ratios as function of s for (a) B — > Kfi'^fi~ 
and (b) B Kt^t~ . The curves with and without resonant shapes represent 
including and no LD contributions, respectively. Legend is the same as Figure 
2. 

Figure 4: The differential decay braching ratios as function of s for (a) B K*fi^fi^ and 
(b) B K*T^T~ . Legend is the same as Figure 3. 

Figure 5: Forward-backward asymmetries for (a) B K*fi^fi~ and (b) B K*t^t~ . 
Legend is the same as Figure 3. 

Figure 6: Longitudinal polarization asymmetries for (a) B K^^^~ and (b) B 
Kt^t~ . Legend is the same as Figure 3. 

Figure 7: Longitudinal polarization asymmetries (a) B K*^^fi^ and (b) B — » K*t^t^ . 
Legend is the same as Figure 3. 

Figure 8: The longitudinal polarization asymmetry for B — > K*t^t^ in the PQCD (solid 
curves) and LF (dashed curves), respectively. 
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Figure 1: Form factors of Fi (solid curves), F2 (dashed curves) and Ft (dotted curves) for 
the B ^ K transition as a function in the PQCD (bold hues) and QM (un-bold hues), 
respectively. 
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Figure 2: Form factors of Fi, F2 and Ft for the B —>■ K transition as a function with 
the dependent wave function in Eq. (|46| ) in the PQCD (sohd curves), QM (dot-dashed 
curves), and LCSR (dashed curves), respectively. 
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Figure 3: The differential decay bracfiing ratios as function of s for (a) B — > Kix^ ix~ and 
(b) B Kt^t~ . The curves with and without resonant shapes represent including and no 
LD contributions, respectively. Legend is the same as Figure 2. 
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Figure 4: The differential decay braching ratios as function of s for (a) B — > K*iji'^iJ, and 
(b) B —>■ K*T^T~ . Legend is the same as Figure 3. 
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Figure 8: The longitudinal polarization asymmetry for B — > K*t^t in the PQCD (solid 
curves) and LF (dashed curves), respectively. 
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